Theory of the quantum Hall effect in graphene 
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We study the quantum Hall effect (QHE) in graphene based on the current injection model. In our 
model, the presence of disorder, the edge-state picture, extended states and localized states, which 
are believed to be indispensable ingredients in describing the QHE, do not play an important role. 
Instead the boundary conditions during the injection into the graphene sheet, which are enforced 
by the presence of the Ohmic contacts, determine the current-voltage characteristics. 

PACS numbers: 73.63.-b,73.21.-b,73.43.-f 



Several experiments have studied the quantum Hall 
effect (QHE) in graphene Recently the observa- 

tion of the fractional QHE in a suspended high-mobility 
graphene device has been reported using a two-terminal 
measurement setup [3, Surprisingly, four-terminal 
measurements in small graphene devices did not reveal 
the QHE. This observation has been related to the large 
influence of the metallic contacts on the formation of the 
Hall potential in small graphene flakes, where the cur- 
rent and voltage probes short out the Hall device 
For the theory of the QHE, the experiments provoke fun- 
damental questions about the role of the metallic con- 
tacts, the correct incorporation of boundary conditions, 
and the importance of electron-electron interactions in 
graphene. These questions are not addressed in the com- 
monly used edge-state and disorder models of the QHE, 
which have been proposed for the QHE in graphene and 
have been taken over from models for conventional semi- 
conductors ^6]. Here, we give further theoretical consid- 
erations which stress the importance of boundary effects 
in small Hall devices and the role of electron-electron in- 
teractions for the integer and fractional QHE. 
Before we discuss the boundary effects, let us see why 
boundary effects are in general not considered to be part 
of theories of the QHE. Two theoretical models are com- 
monly invoked to explain the IQHE: the disorder model 
and the edge-state picture [zl . The disorder model adapts 
a specific distribution of the density of states (DOS) of 
a quasi two-dimensional electron gas in the presence of 
a strong perpendicular magnetic field and a randomly 
fluctuating potential. It is assumed that the DOS splits 
into two parts centered around each Landau level: an 
extended state band at the center of each Landau level 
which is bordered by a broad region of localized states. 
These results are obtained by using ensemble averaged 
Green functions representing an infinite two-dimensional 
system Q. The QHE is viewed as a phase transition, 
which is not related to the measurement apparatus with 
its contacts. The second model of the QHE, the edge- 



state approach, restricts the two-dimensional plane by 
two edges between semi-infinite leads. The strong mag- 
netic field results in a quasi one-dimensional transport 
along the two edges and thus gives rise to two oppo- 
sitely flowing currents. A four terminal measurement 
should guarantee a clean signature of the QHE even in 
the presence of disorder within the device [9| . Both mod- 
els use contrived translational invariances which are not 
present in an actual Hall device. The metallic nature of 
the Ohmic contacts at the device border and the speciflc 
sample geometry play no special role and are ignored. 
This is not the case for the classical Hall effect, where 
the boundary conditions and the device geometry are 
crucial for the calculation of the self-consistent Hall po- 
tential U , 12] . For the determination of the classical 
Hall potential, the metallic contacts have to be consid- 
ered as equipotential surfaces, which enforce also in the 
two-dimensional subsystem a uniform potential under- 
neath the contacts [13]. The two-terminal resistance of 
a classical Hall device can be readily calculated by the 
ratio of the source-drain voltage to the source-drain cur- 
rent. The specific device geometry and the placement 
of metallic contacts do have a strong infiuence on the 




FIG. 1: Hall potential (color online) obtained by solving the 
Laplace equation under the Hall boundary conditions of a 
steady current flow in the presence of a strong magnetic fleld. 
At the upper left corner electrons enter the device in a region 
of high electric fields and move to the lower right corner. 
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Hall potential solution [l3[ . The fundamental reason for 
the formation of the Hall potential are the interactions 
between the electrons in the complete device, including 
the contacts, which translate the magnetic Lorentz force 
acting on every electron into a global, non-trivial adjust- 
ment of the potential with the emergence of two hot-spots 
at opposite corners of the device. The interactions are 
not present in a Fermi liquid model of effectively non- 
interacting electrons [1] and thus these models are not 
sufficient to explain the experimental observations of hot- 
spots and the emergence of the classical Hall field in the 
QHE. 

In the following we describe the injection model of the 
QHE in graphene for a two-terminal measurement. We 
incorporate the metallic contacts in the model and thus 
employ different boundary conditions than either the 
edge-state or disorder models. Our perspective is analo- 
gous to the theory of scanning tunneling microscopy and 
quantum point contacts, in which electron flow in a re- 
stricted region ultimately determines the conductance. 
The mean field potential (Fig. [1]) gives rise to an injec- 
tion hotspot and an exit hotspot, where the drift velocity 
takes on its larges value due to the fast change in the Hall 
potential over a small corner region. Interestingly, images 
of the Hall potential in GaAs/AlGaAs heterostructures 
in Ref. closely resemble Fig. [1] The rate of electrons 
entering the device within the hotspot region is given (as 
in the edge-state model) by the convolution integral of 
the LDOS at the injection points with the group velocity 
and the Fermi-Dirac distribution. The main difference 
to an edge-state model is the absence of translational 
invariance due to the presence of the contacts and the 
strong bias of the current injection towards a corner of 
the device. The resulting current fiows unidirectional di- 
agonally across the device. 

From the classical solution alone, no reason for the quan- 
tization of the resistivity exists. The quantization re- 
quires to study the relation between the source-drain 
voltage and the current in a Hall bar. The self-consistent 
potential in a Hall device does depend on the applied 
voltages, the magnetic field, and the geometry of the 
sample. For rectangular samples with a length to width 
aspect ratio L/W ^ 1, the solution of the Laplace equa- 
tion for a Hall bar in the presence of a magnetic field and 
current leads to a geometry independent solution with 
high electric field in two opposite corners of the device 
[ill] , where the Hall potential attains the universal form 



lution as the mean-field potential, which emerges due to 
interaction and screening between the electrons and the 
positive charges and by considering the metallic bound- 
ary conditions at the current source and sinks [l3|. We 
study the propagation of the effectively non-interacting 
electrons in the mean-field potential. We view our model 
as providing a starting point for a more rigorous inclusion 
of interaction effects, already in the non-fractional quan- 
tum Hall effect. Experimental evidence for the existence 
and relevance of the mean-field solution is provided by 
the absence of the QHE in a four-terminal measurement 
in graphene, which has been attributed to the dominance 
of the metallic boundary conditions at the contacts in 
small devices j^]. For low energies, the nearest- neighbor 
tight-binding Hamiltonian of graphene can be expressed 
by an effective Dirac-type Hamiltonian . Including the 
minimal coupling of the potential V{x, y) and the mag- 
netic field via the vector potential A. = B{—y, 0, 0) to the 
kinematic momentum IT = p — yields: 
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The four components of the wavefunction are labelled 

(V'i,i/'2,'03,''/'4) = (V'a >V'f >A">V'f' ), where K and K' 
refer to the two K points in the first Brillouin zone and 
the two sublattices A and B, and c w 10^ m/s denotes 
a velocity. The Hamiltonian omits the spin degree of 
freedom, which is added later using an effective g-factor. 
The propagation of wave packets by numerical methods 
allows us to accurately determine the local density of 
states (LDOS) without the need to introduce half-infinite 
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Here, Vsd denotes the voltage difference between the 
source and drain contacts. For the inverse geometrical ra- 
tio, L/W ^ 1 (used for measuring the FQHE in graphene 
[3|), a very similar solution emerges, since the long con- 
tact region enforce the boundary conditions very effi- 
ciently. In the following, we will use the uniform-map so- 



FIG. 2: Overview of the LDOS (color online) with an effec- 
tive jr-factor g* — 2 for the electric field £ = 200 kV/m and 
the magnetic field Z? = 15 T as a function of the energy in 
Eq. (|13p . The red and blue areas are the contributions from 
the spin-up and spin-down components. Note the modulation 
within each spin-split component, which is caused by the elec- 
tric field and increases for higher levels. The black solid line 
shows pxy = The energy is given in units of 141 meV. 
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leads, which are not compatible with the boundary con- 
ditions discussed above. We have adapted the approach 
of Ref. [l^ to the massless Dirac equation in graphcnc 
using a recursively evaluated polynomial expansion of the 
time evolution operator [l^. We calculate the LDOS in 
a strong magnetic field and for the potential given in 
Eq. IT]) by tracking the time-dependent autocorrelation 
function for several picoseconds. The resulting LDOS is 
intrinsically broadened (by decay in the autocorrelation 
function due to flux leaving the injection area under the 
combined influence of the electric field potential and the 
magnetic field) and shows only extended states which 
connect one corner of the device with the opposite one. 
We find that our numerical results for the LDOS in the 
corner potential ^ are in excellent agreement with the 
analytically derived LDOS in perpendicular and homo- 
geneous electric and magnetic fields, provided we choose 
the homogeneous electric field value to match the local 
potential gradient: 
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The uniform field case is analytically solvable using the 
proper-time approach of Fock [l^, El, El, HO] . We con- 
struct the LDOS in the magnetic field B = (0, 0, B) and 
the electric field £ = (0,£,0) from the four scalar com- 
ponents and the eigenenergies 
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given in Landau gauge [2l|. The spinor eigenfunctions 
can be written in terms of oscillator functions 
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Using the coefficients 
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with the plane wave solution in the a;-direction, (/)p^ (x) 
we obtain 
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ipi,n,pAr) = 4>pjx)[ —j (sgn(n) aM|„|(0-bw[„|_i(^)), 

(8) 



■03,n,p,(r) = V'2,n,p^(r), V'4,n,p, (r) = V'l^n.p.lr). (10) 

For n — 0, the M|„|_i(^) terms vanish and the wavefunc- 
tions acquire an additional normalization factor of ^/2. 
We obtain the LDOS by evaluating 
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In the presence of an electric field, the spinor index i 
can no longer be used to identify the sublattice. The 
introduction of the real spin completes the computation 
of the LDOS 
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where p,B denotes the Bohr magneton (which contains 
the normal electron mass, not the effective one) and g* is 
the effectice g-factor of the electron. The LDOS shown 
in Fig. [2] is symmetric with respect to the E — Q value. 
The centers of the spin-split LDOS are located at en- 
ergies En^p=o ± ^g*^B- For electric fields approaching 
£/B — c, the harmonic oscillator functions are replaced 
by Airy functions [19,]. The transition from a magnetic 
field dominated LDOS to an electric field dominated one 
for the non-relativistic case is also discussed in Ref. [2^ . 
The group velocity dEn^p^/dpx in the relativistic case is 
given by £/B and thus unchanged in the non-relativistic 
limit and does not depend on the energy of the particle. 
In the injection model, the LDOS in the injection region 
together with the group velocity determines the injected 
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FIG. 3: Hall resistivity as function of Fermi energy (color 
,e/3\ online). Shown are the filling factors = 2 — 10 at a tern- 

V'2,«,pjr) = (f>p^{x) ( — j (a M|„|_i(e)-sgn(n) b ii|„|(0), perature of T = 5 K for two dc-currents Isd = 1 M and 

(■9) Isd = 10 fiA. 
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current density. The magnitude of the electric field in 
the middle of the device is given by Vsd/W and thus a 
higher source-drain voltages increases the drift velocity. 
An important quantity is the statistical distribution of 
the electric field strengths at the injection sites. Since 
fully quantum-mechanical simulations of open quantum- 
systems with Coulomb interactions are not available, we 
have to make assumptions about the magnitude and spa- 
tial distribution of the electric field values. Also dissipa- 
tive processes are expected to occur at the hot-spots at 
the two opposite corners with the highest probability for 
electrons to enter and to leave the graphene flake. The 
superposition of various field values leads to a suppres- 
sion of the electric field induced gaps within the Landau 
levels. The current is given by the product of the group 
velocity and the LDOS convoluted with the temperature- 
dependent Fermi-Dirac distribution /: 

IsD ^eJ^dEj dre f{E - Ep)^nl^^-'' {r,; E), 

(14) 

where the magnitude of the electric field in the injec- 
tion region near the hot-spot depends on the value of 
the source-drain voltage drop occurring in this corner, 
given by Eq. ([3]). Thus we obtain a voltage dependent 
group velocity and broadening of the LDOS, depending 
on the injection points at Tc- For a theoretical simulation 
of the QHE in constant current mode, we have to iter- 
ate Eq. (|14p with different values of Vsd until we match 
the desired total current. In Fig. [3] we show the cur- 
rent dependence of the Hall curves for a magnetic field 
of ;B = 10 T as a function of Fermi energy, which can be 
adjusted experimentally by varying the voltage of a back- 
gate. The figures are calculated under the assumption, 
that a single electric field strength near the hot-spot dom- 
inates the current, which is 50 times stronger than the 
linear Hall field in the middle of the device. The increase 
of the local-electric field at the injection region of a fac- 
tor 50 is an estimate based on the extension of observed 
hot-spots in GaAs/AlGaAs heterostructues For cur- 
rents in the order of 10 /iA a strong modulation of the 
Landau levels is visible which develop a very asymmetric 
shape and a gap near the filling factors v = 4,8, 12, . . ., 
where the resistivity curves intersect at current indepen- 
dent points between two adjacent plateaus. These modu- 
lations and gaps are not caused by interactions but by the 
peculiar shape of the LDOS in graphene in the current 
injecting corner. We expect that these hot-spot induced 
gaps are reduced by the simultaneous emission from the 
other injection points with lower electric field values. 
In conclusion, we have calculated the source-drain- 
voltage-current relation of a graphene device in strong 
magnetic fields. We have developed a theory of the QHE 
in graphene, which shows a rich substructure and modu- 
lation of Landau levels as a function of the back-gate volt- 
age due to the presence of hot-spots in the device. The 



modulation in graphene is different compared to other 
semiconductors due to the perfect alignment of the nth 
and {n — l)th Landau level around the same energy. The 
experimental observation of current induced structures 
can clarify the electronic transport paths in graphene and 
pin down the region where electrons enter a device in the 
presence of high magnetic fields. The presence of disor- 
der is not a requirement for the existence of the QHE. 
Puddles of localized electrons away from the current in- 
jecting corner are not able to stop the injection process 
and thus may only deform the pathways of the electrons 
through the sample, which affect the longitudinal voltage 
drop, but not directly the Hall conductivity. 
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of the DFG (KR 2889-2/1). We thank P. Kramer, 
C. F. Huang, L.-H. Lin, J. Fabian, and D. Weiss for help- 
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